Possible vortex splitting in the cuprate superconductors 
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We propose that the observed sphtting of the vortices in the cuprates into fractional vortices 
(partons) may be of static rather than of dynamic origin. This interpretation is backed by a study 
of a model with a dominant d-wave and subdominant s-wave pairing interaction. We find that the 
vortex may split into two partons, both of which carry one half of the magnetic flux quantum. The 
partons are hold together by a confining string along which the phase jumps approximately by n 
and their equilibrium distance increases with lowering the energy difference e between the pairing 
states. The partons become deconfined at the critical point where e vanishes. 



The nonsuperconducting phase of the high temper- 
ature superconductors exhibits anomalous features 
On the other hand, the low-temperature superconduct- 
ing state is believed to be well described by the standard 
Bardeen-Cooper-SchriefFer paradigm, if the d-wave pair- 
ing symmetry and the Landau Fermi liquid corrections 
are taken into account. In view of the anomalies of the 
normal phase, it is tempting to look for unconventional 
features in the superconducting phase as well. Interesting 
results have in fact been obtained by the scanning tun- 
neling spectroscopy of the vortices in Bi2Sr2CaCu208+(5 
@) S] J according to which the vortex cores may split into 
several subcomponents with a spacing in the range 10 
- 100 A. The experimental results were interpreted as a 
result of the vortex hopping between different pinning 
sites 0, 01 • The purpose of this paper is to propose an 
alternative interpretation of the experimental results, in 
which the splitting of the vortices into partons is consid- 
ered to be of static origin. Based on an analogy with the 
physics of dislocations, we will show that the vortex may 
split into two fractional vortices (partons), each of them 
carrying one half of the magnetic flux quantum $0- 

It is well known that screw dislocations in fee materi- 
als can split into two Shockley partial dislocations whose 
(fractional) Burgers vectors add up to an integer lattice 
vector [5| . In that case the singular dislocation line trans- 
forms into a singular strip whose borders are formed by 
the partial dislocations. The relative displacement of the 
crystal on both sides of the strip is not equal to a lattice 
vector. The dynamical reasons for the stability of the 
partial dislocations are: (i) repulsion between the paral- 
lel partial dislocations and (ii) the low elastic energy cost 
of the displacement across the strip. In the vortex case 
we will show that in addition to (i), which is always true, 
the criterion (ii) may be satisfied in superconductors with 
sufficiently strong subleading pairing interactions. 

Besides serving as an alternative explanation of the 
experiments 0, the parton hypothesis provides addi- 
tional support to the interpretations of the pseudogap in 
the high-temperature superconductors as an incoherent 
liquid of singlet electron pairs on the bonds of the Cu02 
lattice The major open problem in this line of think- 



ing is the question about the mechanism leading to the 
phase disordering of the pairs. It has been argued that in 
order to destroy the phase ordering and to stabilize the 
pseudogap state, the presence of cheap vortices in the 
cuprates is required [l|]. We will show that the energy 
of the split vortices may be substantially reduced with 
respect to the conventional vortex energy. 

Our main assumption is the existence, in addition to 
the leading d-wave interactions, of subleading pairing in- 
teractions in the s-wave sector. Our motivation is as 
follows. It seems reasonable to assume that the model 
of the cuprates should contain a strong on-site repulsion 
and a moderate antiferromagnetic nearest neighbor spin- 
spin interaction. It is well known that within this type 
of a model, condensates with both cZ-wave and s-wave 
symmetry may form Q. On the other hand, we are not 
aware of direct experimental evidence for such subleading 
pairing tendencies in the cuprates. It has been argued, 
however, that the large second harmonics of the current- 
phase relation observed in the cuprate grain boundary 
Josephson junctions [1, Q provides an indirect evidence 
for the existence of subleading pairing interactions [10]. 

Previous works have noted that the structure of iso- 
lated vortex lines in superconductors with competing 
pairing interactions may become very rich. In particular, 
within phenomenological Ginzburg-Landau theory it has 
been shown that in the vicinity of the vortex cores in d- 
wave superconductors, there may nucleate a finite s-wave 
component with a nontrivial phase structure fll'|. How- 
ever, these results do not explain the experiment since 
the dominant d-wave field has only one singularity. For 
the same reason, neither the more recently proposed non- 
singular vortices [l^ can explain the experimental data. 
Also Volovik has suggested [13| that vortex splitting may 
have been observed in [5], but he has not presented any 
calculation to support this hypothesis. 

We model the Cu02 plane as a square array of super- 
conducting islands. There are two alternative interpreta- 
tions of this lattice. Phenomenologically, one may think 
of it as a coarse-grained model of the Cu02 plane. Micro- 
scopically, adopting the short-range RVB picture for the 
sake of simplicity, the lattice may be thought of as the 
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set of the centers of mass of the nearest-neighbor Cu-Cu 
singlets. The fluctuations of the superconducting ampli- 
tudes of the islands are neglected and the only dynamical 
variable describing the island at i is supposed to be the 
phase of the condensate 6i. The islands are assumed to 
be coupled by the Josephson effect and we postulate that 
the Hamiltonian of the plane reads as 

H [^cos{e, - Oj) - Jcos(26'i - 29^)] , (1) 

where the sum is taken over the nearest-neighbor sites. 
The d-wave pairing state is described by £ > and in 
the spin language it corresponds to an antiferromagnetic 
configuration of 9i. Note that for e < it is the s-wave 
pairing state that is stable and therefore e = corre- 
sponds to a quantum critical point. We assume that 
J > and therefore the ground state phase difference 
in Eq. ((T|) jumps discontinously from tt at e > to at 
£ < 0. The alternative choice J < would correspond 
to a continuous change of the ground-state phase differ- 
ence, thus physically corresponding to a d+is state in the 
vicinity of £ = 0, i.e. to a time reversal-breaking mixture 
of the d-wave and s-wave pairing states. Such cooper- 
ation of different pairing states is generically favourable 
at weak couphng [3]. On the other hand, our case J > 
corresponds to a competition between the pairing states. 

In order to proceed we modify the model Eq. ([1]) in 
several ways. First, we perform a gauge transformation 
&i — > 9i + on one of the sublattices. This changes 
the sign of the first term in Eq. ^ and, as a result, 
the d-wave state corresponds to the ferromagnetic state. 
Second, we redefine the zero of energy so that the ho- 
mogeneous case corresponds to E = and finally, we 
include the coupling to the magnetic field. For the sake 
of simplicity, we consider a layered tetragonal material 
with in-plane lattice constant d and c-axis lattice con- 
stant dc and we consider only vortices along the c axis. 
Finally we have 

E-j:^{0.,) + l^Y.^l (2) 

where E — E / (AJ + e) is the dimensionless vortex energy 
per Cu02 plane and A"^ = in^ fioi'^J + £)/(*o^c)- 

The first term in Eq. ^ corresponds to the sum of 
dimensionless Josephson energies of single bonds, 

e{9) = (1 - 2c)(l - cos6l) + |(1 - cos 26*), (3) 

where c = 2J/(4J -I- £) is a parameter measuring the 
strength of the second harmonic contribution to the 
Josephson energy of the bonds. We have introduced a 
dimensionless vector potential atj — 2t: f!' A ■ dr/$o and 
a gauge invariant phase difference between lattice sites i 
and j , 6ij = 9j — 0i + aij . 



The second term in Eq. ^ corresponds to the energy of 
the magnetic field, ipt is the dimensionless flux threading 
the plaquette with lower left point at i. If the plaquette is 
formed by the points ijkl, then Lpi — aij + ajk + aki + an. 

Minimizing the energy Eq. ([2|) with respect to we 
obtain the coupled set of discretized Maxwell equations 

ift-y-iPi = ■J^j{9.i + x - 9t + Gi^t + x), (4) 

d^ 

where x and y are elementary lattice vectors in the x 
and y directions, respectively, and we have introduced 
the dimensionless current j{9) = (1 — 2c) sm9 + csin26'. 
Note that for slowly varying fields, 16*1 <C f, = 6* for 
all values of c. Making use of Eqs. (|4l5p . we thus identify 
A as the penetration depth and from A « 2600 A [l3| and 
dc ~ 7.5 A, we estimate 4 J -|- £ « 6 meV independently 
of c. Thus c is the only free parameter in the theory. 

In what follows we make use of the inequality A S> c? 
and solve Eqs. (|4|5p perturbatively with the small param- 
eter d/X. To this end let us require that the phase field 
9i satisfies the following equations for all sites i: 

^j(^.+r-e.)=0, (6) 

T 

where the sum is taken over the four nearest neighbor 
directions r. Physically this corresponds to a lattice ver- 
sion of the continuity equation V • j = and more for- 
mally it might be thought of as a discretized version of 
the equation = 0. Once Eq. ([6]) is satisfied, one can 
estimate the magnetic fields by neglecting the vector po- 
tential Gij on the right-hand sides of Eqs. (|4I5|) and one 
finds that (fi oc d^/A^ <C 1. The vector potential at dis- 
tance R from the vortex center can therefore be chosen as 
a ~ {d/X)^{R/d). On the other hand, the typical phase 
difference 9i+T — 6i at distance R is d/R. Therefore the 
current distribution is well described by the phase-only 
solution inside the circle with radius R X. Beyond this 
range the correct solution should differ only marginally 
from the standard vortex solution [l^. In this paper we 
therefore concentrate only on the region R <ti X. 

Our task is therefore to find a solution to Eq. ([6]) with a 
finite winding number. Consider first the standard solu- 
tion describing a phase winding by 2tt around the vortex 
center at (0.5,0.5). Consider further the straight fine 
passing through the vortex center and parallel to the x 
axis. This line cuts the set of bonds connecting points 
{x, 0) and {x, 1) of the lattice. We will work in a gauge 
where the phase difference A9x ~ 9xo — 9xi changes be- 
tween 2tt for X — > — oo and for x +oo. 

We have solved Eq. ([6]) on lattices 200 x 200 numerically 
by the standard iterative procedure. By varying the ini- 
tial configuration, several solutions could be found. The 
number of different solutions increased with c. For all 
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FIG. 1: Phase difference /S.9x ~ 9xO ~ d^i for bonds which 
are cut by a straight line passing through the vortex center 
and parallel to the x axis, plotted versus the position x of 
the bond. The partons are located at a:: = 0.5 ± a. Note the 
growth of the plateau at A9x = tt as the quantum critical 
point c = 0.5 is approached. 



solutions we calculated their Joscphson energy and for 
every studied c, we have identified the optimal solution 
with minimal energy. From now on, we will focus on the 
optimal solutions. 

The results for the phase jump AO^ for several values 
of c are shown in Fig. [TJ In agreement with our expecta- 
tions, as a function of x, Ad^ exhibits only a single step 
[li^ for sufficiently small c. However, for c > 0.25 two 
partial steps develop. The phase jump AO^ at those steps 
changes from to vr and from tt to 2tt, respectively, cor- 
responding to the presence of two partons. The spacing 
between the partons, 2a, defined as the distance between 
the steps, is plotted in Fig. [2] Note the steep increase 
of a(c) for c approaching the quantum critical point at 
c = 0.5. 




FIG. 2: The equilibrium half-distance a between the partons 
(in units of d) as a function of c. The inset shows the energy 
of all solutions found for c = 0.475 as a function of their a. 
On the vertical scale, only energy differences are meaningful, 
since the energy contains an additive constant dependent on 
the lattice size. 

Once the phase fields are known, the magnetic field 
distribution can be calculated making use of Eqs. (I4I5I) . 
The result of this calculation is shown in Fig. [3] which 
explicitly shows that the vortex splits into two partons. 



FIG. 3; Spatial distribution of the dimensionless magnetic 
flux ifi in a split vortex for c = 0.475. 
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FIG. 4: Kinetic energy (in units of ij + e) of the bonds along 
the X direction (a) and along the y direction (b) for c = 0.475. 



Let us discuss the energetics of vortex splitting in more 
detail. Figure |4^ shows the Josephson (kinetic) energy of 
the bonds in the x direction. By the standard argument 
about instability of multiply charged vortices one can see 
that this configuration has a lower total kinetic energy in 
the X direction, E^, than the usual vortex solution. Thus, 
the gradient of Ex pushes the partons apart [l^ • 

There is however an opposing force, which is caused 
by the necessary existence of a cut joining the partons, 
across which the phase has to jump approximately by 
TT. The kinetic energy of the bonds in the y direction 
is shown in Fig. |3|3. One can see explicitly that, at long 
distances 2a between the partons, the total kinetic energy 
in the y direction, Ey, grows linearly with a due to the 
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energy cost of the cut, thus forming a confining potential 
for the partons. 

The competition between the repulsive and attractive 
forces leads to the presence of a minimum of the total 
vortex energy as a function of the parton distance, as 
demonstrated explicitly in the inset to Fig. O 

So why don't the partons always form? The reason for 
this is that, for c < 0.25, the Josephson energy Eq. ([3]) 
has a local maximum at = tt. This destabilizes the 
parton solution at those values of c. For c > 0.25, the 
parton is locally stable, but the equilibrium interparton 
distance is in general very short, see Fig. [2l because the 
string tension associated with the cut is huge. However, 
when approaching the quantum critical point c = 0.5, the 
energy difference e(7r) — e(0) measuring the string tension 
diminishes and the interparton distance 2a grows. 

The theory as developed so far applies only to situa- 
tions when a ^ A. Therefore it does not apply in the 
immediate vicinity of the quantum critical point. In or- 
der to check the robustness of our picture we have solved 
the full equations ()4I5|) at the critical point c = 0.5 for 
a single fractional vortex. To this end we have rewritten 
Eqs. (|4l5)) as 



2X\ 



+ 7r5,o, (7) 



where the first term can be thought of as a discretized 
Laplacian and only the nearest-neighbor sites contribute 
to the sum over r. The second term is the source term for 
a vortex whose phase winds by tt. The solution to Eq. ([7|) 
was obtained by the standard iterative procedure. We 
have checked explicitly that in the core region i? ^ A the 
approximate solution which satisfies Eq. ([5]) is in perfect 
agreement with the full solution. 

The conventional vortex energy at c = can be 
estimated using the London theory [l3| and we find 
i?(0) « 7rln(A/(i). Taking d in the interval between 
3.8A/\/2 and C w 14A [18|, we estimate £'(0) w 16-22. 
From our numerical solution we can determine the change 
of the vortex energy with c and, e.g. at c = 0.475, 
^(0) - i?(0.475) w 4.85. It is seen that finite values of c 
may lead to a substantial reduction of the vortex energy. 

Before concluding let us discuss the relevance of our 
results to the experiment 0. Following Q, we assume 
that the pinning effects are decisive in determining the 
vortex shape. The difference with respect to Q is that in 
our picture all vortices in a perfect sample should be split. 
The observation of unsplit vortices can be explained by 
the attraction of the partons to the same pinning center 
by a potential which is stronger than the energy gain 
due to splitting Moreover, the interpretation of the 
observed slow temporal evolution of the vortex shapes 
proposed in Q is applicable also in our picture. 

In 0], the quantum tunneling of vortices between 



nearby pinning sites was considered as the most likely 
mechanism of vortex delocalization. Thus the vortex 
should be described by a linear superposition of wave- 
functions describing the vortex localized at the various 
pinning sites. It is a subtle issue to distinguish a vortex 
described by such a wavefunction from our stable parton 
picture. We believe the best way to distinguish these two 
alternatives is to determine whether split vortices form 
also in perfect samples. 

In conclusion, within a simple model we have shown 
that, in the vicinity of a quantum phase transition be- 
tween pairing states of different symmetry, the vortex 
cores may acquire a nontrivial parton structure. We be- 
lieve that this result is interesting in several respects: (i) 
in scanning tunneling spectroscopy of the cuprates, the 
parton structure of the vortex cores has been observed 
in some of the samples; (ii) the phenomenon is analogous 
to the formation of the partial Shockley dislocations; (iii) 
the parton structure lowers the energy of the vortices and 
this may be relevant for the interpretation of the pseu- 
dogap phase; (iv) precisely at the quantum critical point 
e = 0, fractional vortices become deconfined. 
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